Charge-density waves in the Hubbard chain: evidence for 4k p instability 
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Charge density waves in the Hubbard chain are studied by means of finite-temperature Quantum 
Monte Carlo simulations and Lanczos diagonalizations for the ground state. We present results 
both for the charge susceptibilities and for the charge structure factor at densities p = 1/6 and 1/3; 
for p = 1/2 (quarter filled) we only present results for the charge structure factor. The data are 
consistent with attp instability dominating over the 2k,F one, at least for sufficiently large values 
of the Coulomb repulsion, U. This can only be reconciled with the Luttinger liquid analyses if the 
amplitude of the 2&f contribution vanishes above some U*(p). 
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•4-* 



C 

O 

o 



m 
O 

(N 
O 
O 

o 

-I— > 



T3 

c 

o 
o 



x 



Charge-density waves (CDW's) are present in a vari- 
ety of strongly correlated electron systems, ganging from 
quasi-one-dimensional organic canductorsErB to the more 
recently discovered manganites.Q In order to understand 
the influence of CDW formation on magnetic order and 
transport properties, a crucial issue is to establish the 
period of charge modulation in the ground state. This 
period, in turn, is primarily determined by the inter- 
play between electron-phonon and electron-electron cou- 
plings. For the specific example of quasi-one dimen- 
sional organic conductors at quarter-filled band, both, 
period-2 and period-4 modulations have been observedp 
these correspond, respectively, to Akp and 2kp , where 
kp = up/2 is the Fermi wave vector for a density p of 
free electrons on a periodic lattice. The use of simpli- 
fied effective models capturing the basic physical ingre- 
dients should therefore be extremely helpful in predict- 
ing the dominant instability.!] In this context, the Hub- 
bard model can be thought of as a limiting case (of van- 
ishing electron-phonon interaction), in which the influ- 
ence of electronic correlations on CDW modulation can 
be monitored. However, even for this simplest possible 
model there has been a disagreement between analyses 
of the continuum (Luttinger liquid) version,ETu and early 
finite temperature .(world-line) quantum Monte Carlo 
(QMC) simulations.L4J According to the former, the large- 
distance behaviour of the charge density correlation func- 
tion is given byl3 
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where the amplitudes A\ and A2, and the exponent K p 
are interaction- and density-dependent parameters; for 
repulsive interactions!] | < K p < 1, so that charge cor- 
relations are expected to be dominated by the 2kp term. 
By contrast, the simulations pointed towards Akp being 
the main correlations. Nonetheless, based on-a Renor- 
malization Group (RG) analysis, it was arguedt2l that the 



2kF instability should eventually dominate over the Akp 
one for sufficiently low temperatures. For infinite cou- 
pling the system becomes effegt^ely a spinless fermion 
problem with a Akp instabilityJZrEdO and there is no dis- 
agreement in this case. 

Since present day computational capabilities allow one 
to reach much lower temperatures and larger system sizes 
than before, a numerical reanalysis of the model is cer- 
tainly in order. Our purpose here is to present the results 
of such a study. The Hubbard Hamiltonian reads 
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where, in standard notation, U is the on-site Coulomb 
repulsion; the hopping integral sets the energy scale, so 
we take t = 1 throughout this paper. We probe fi- 
nite tempeijatuEe properties through determinantal QMC 
simulationsEaila for the grand-canonical version of (0), 
H = Ti — fiN, where N = ^Zia n ^ the chemical po- 
tential pi is adjusted to yield the desired particle density. 
This analysis is supplemented by zero temperature calcu- 
lations: the ground state of Eq. (||), for finite lattices of 
N a sites with periodic boundary-conditions, is obtained 
through the Lanczos algorithm,E3ila in the subspace of 
fixed particle-density (canonical ensemble). 

The signature of a CDW instability is a peak at q = 
q* in the zero-temperature limit of the charge-density 
susceptibility, 
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where the imaginary-time dependence of the operators is 
given by rii(r) = e rW n,; e~ rW , with — + n^. We 
recall that in simulations 'time' is discretized in intervals 
Ar, such that the size along this direction is L = /3/Ar. 
The CDW instability should also show up as a cusp, again 
at q = q* , in the zero-temperature charge-density struc- 
ture factor, 
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where |0) is the ground state. 
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FIG. 1. QMC results for the charge susceptibility as a func- 
tion of the inverse temperature /3 for a chain with N s — 36 
sites and occupation p = 1/6: (a) U = 2, (b) E/ = 6, and (c) 
U = 9. Circles and squares represent data for q — 2kF and 
respectively. Error bars are smaller than data points, 
and all lines are guides to the eye only. 




FIG. 2. QMC results for the charge susceptibility as a func- 
tion of In 7V S for U = 9; each point is obtained at j3 = N a /4. 
Circles and squares represent data for q = 2/cf and 4/cf, re- 
spectively. The straight line through data for 4/cf is a least 
squares fit. 



Let us first discuss results for an electronic density 
p = 1/6, for which we found that the fermionic determi- 
nants in QMC simulations do not suffer from the 'minus- 



sign problem' this allowed us to reach inverse tem- 
peratures as large as (3 — 25. ^From Fig. |l| we see that 
the 4/cf charge susceptibility appears to be increasing 
with decreasing temperature, at a rate faster than that 
at 2k f, especially for U — 6 and 9. The data in Fig. 
|l| were obtained for At = 0.125, but we have explicitly 
tested other values to ensure they do not change signif- 
icantly as At — * 0; each datum point involves typically 
20,000 QMC sweeps over all time slices. Further, in order 
to check if this increase is limited by finite-size or finite- 
temperature effects, we performcd-additional simulations 
on a 'square space-time lattice' ,E9 i.e., we set N s = L; 
the result is displayed in Fig. |^ for N s = L < 96. While 
the charge susceptibility at 2kp seems to saturate as N s 
increases, the one at 4/cf still scales with ln7V s , up to 
the largest sizes considered. Thus, in spite of the very 
low temperatures reached, we were still unable to find 
indications of a crossover temperature below which the 
system is dominated by the 2kp instability. Also, since 
N s oc l/T for the data in Fig. ||, the 4fcp charge suscepti- 
bility grows logarithmically with the temperature in this 
range, similarly to the infinite coupling limit.E3 
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FIG. 3. Charge structure factor for electron density 
p — 1/6 and U — (up triangles), 2 (left triangles), 4 
(squares), 6 (down triangles), 9 (diamonds), 12 (circles), and 
20 (right triangles). The system size is N s = 24, and succes- 
sive vertical shifts by -0.02 have been imposed on the curves, 
for clarity. 



We now discuss the charge structure factor at zero tem- 
perature, as obtained from Lanczos diagonalizations, still 
for p = 1/6. Figure || shows C(q) for several values of 
U; for clarity, the curves are shown after suffering suc- 
cessive displacements. For the free case, U = 0, we see 
a sharp plateau beginning at q = 2kF = 7r/6, which is 
the signature of the Peierls instability. This behaviour is 
quite different from the one observed for U ^ 0: though 
somewhat rounded for U — 2 and 4, the plateaux now 
start at q = AkF = 7r/3. It is instructive to examine how 
these roundings evolve with system size. In Fig. ^, we 
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single out data for C(q) with [7 = 3 and 12, and sizes 
N s = 12 and 24. For each value of U, the data below 
and above Ak f respectively move down and up as N s 
increases, thus sharpening the cusp; in addition, the po- 
sition of the latter shows no tendency of shifting from 
Akp. Thus, the Lanczos results are consistent with those 
from QMC simulations, in the sense that a Akp instabil- 
ity is dominant, already for moderate values of U. At this 
point, it is worth pointing out that the charge structure 
factor for the Hubbard chain with second neighbour hop- 
ping has been calculated tlppugh density matrix renor- 
malization group (DMRG).t3 Though their interest was 
to extract K p from the slope of C(q) at q = 0, if one 
reinterprets those data along the lines discussed here, a 
predominance of the AkF instability can be clearly in- 
ferred for the largest values of U shown in Fig. 11 of Ref. 
|l9| . DMRG calculations for the two-leg Hubbard ladder 
have also led to 4fc^-like charge correlationsE9 
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FIG. 4. Charge structure factor for electron density 
p — 1/6. Empty squares: U = 3, N B = 12; filled squares: 
U = 3, N s = 24; empty circles: U = 12, N s = 12; filled 
circles: U = 12, N a = 24. A vertical shift by -0.025 has been 
imposed on the curves for U = 12, for clarity. 



We now change the band filling to p = 1/3; the average 
sign of the fermionic determinant is ~ 0.9 in the worst 
cases, thus posing no problems to the resulting averages. 
Figure || shows QMC data for the charge susceptibilities 
with U — 6 and [7 = 8. For U = 6 an upturn at lower 
temperatures seems to be setting in for Akp, while the 
2fcp data show no noticeable change in growth rate. On 
the other hand, for U = 8 the 4kp susceptibility grows 
unequivocally faster with than the one at 2kF', see Fig. 
||. The corresponding Lanczos data for the charge struc- 
ture factor are shown in Fig. ^ and, similarly to Fig. |^, 
the cusp at q = Akp = 27r/3 gets visibly sharper as U 
increases. 




FIG. 5. QMC results for the charge susceptibility as a func- 
tion of the inverse temperature j3 for a chain with N a = 36 
sites and occupation p — 1/3: (a) U = 6 and (b) U — 8. 
Circles and squares represent data for q — 2fc_F and 4&f, re- 
spectively. Error bars are smaller than data points, and all 
lines are guides to the eye only. 
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FIG. 6. Charge structure factor for electron density 
p = 1/3 and U = (up triangles), 4 (squares), 6 (down tri- 
angles), 8 (left triangles), 9 (diamonds), 12 (circles), and 20 
(right triangles). The system size is N s = 18, and successive 
vertical shifts by -0.03 have been imposed on the curves, for 
clarity. 



Unfortunately, for larger band fillings the 'minus-sign 
problem' prevents us from reaching very low tempera- 
tures. Nonetheless, down to the lowest temperatures 
probed with acceptable average signs of the fermionic 
determinant (i.e., (sign) ~ 0.7 at T ~ 1/20), no indica- 
tions of a 2kp peak dominating the charge susceptibility 
were found for p = 1/2 or 3/4. Accordingly, the charge 
structure factor at T = 0, calculated through Lanczos di- 
agonalizations on a 16-site chain at quarter filling, shown 
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in Fig. \A confirms the previous patterns: there is some 
rounding near q — 4kp — tt, which sharpens as U in- 
creases, consistently with a 4kp instability setting in for 
finite f/'s. 



of the amplitude of the 2kp contribution. We hope our 
results stimulate more extensive work both on the Lut- 
tinger liquid and on the lattice model in order to extract 
a quantitative behaviour for the amplitude A\{p, U). 
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FIG. 7. Charge structure factor for electron density 
p = 1/2 and (7 = (up triangles), 3 (squares), 6 (down trian- 
gles), 9 (diamonds), 12 (circles). The system size is N 3 — 16, 
and successive vertical shifts by -0.03 have been imposed on 
the curves, for clarity. 

In summary, for all band fillings examined, the charge 
instability seems to be characterized by a Akp modula- 
tion, rather than by 2kp, at least for U greater than 
some U*(p). The question of how can these findings 
be reconciled with the analyses of the continuum model 
still remains. Since Lanczos data have been obtained 
at zero temperature, and QMC simulations reached 
much lower temperatures than before|L9 the scenario of 
a temperature-driven crossover seems now unlikely; it 
should be recalled that this crossover was predicted based 
on a weak coupling RG analysis. We therefore envisage 
the following scenario. While analyses of the Luttinger 
liquid have so far provided detailed insight into the be- 
haviour of the exponent K p , little is known about the 
dependence of the amplitude A\ of the 2kp contribu- 
tion [see Eq. (]!])] with p and with the coupling constant. 
Our results may be indicating that, for fixed density p, 
A\ — > very fast with increasing U, either exponentially 
or, less likely, as (U*(p)-Uf, for U < U*{p), with 
ip > 1; a crude examination of the roundings near 4k f in 
the structure factors is consistent with U* growing with 
p. An alternative scenario could be that 2kp charge cor- 
relations in the lattice model would suffer from unusu- 
ally slow finite-size effects, thus hindering any present 
day numerical calculations to detect their predominance 
over the 4k f ones; however, since one would need effects 
slower than those suggested by the logarithmic 'correc- 
tion' in Eq. ([!]), this scenario is less appealing. Therefore, 
the presence of a 4k p instability can be made compati- 
ble with Luttinger liquid picture through the behaviour 
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